'''
Created on Aug 9, 2013

@author: csiebern
'''

import numpy as np
import T3D
import matplotlib.pyplot as plt

LASFile = '../LiDAR/Rt_Seiler/Rt_130718(1).las'

LiDAR = T3D.LiDAR(LASFile)

D = LiDAR.data['G3D_DISTANCE']
H = LiDAR.data['G3D_AngleH']*(180/np.pi)
V = LiDAR.data['G3D_AngleV']*(180/np.pi)

plt.figure(1)
plt.subplot(311)
plt.plot(V[0:10000])
plt.title('LAS T3D_AngleV over 10000 points')
plt.ylabel('AngleV [degrees]')

plt.subplot(312)
plt.plot(H[0:10000])
plt.title('LAS T3D_AngleH over 10000 points')
plt.ylabel('AngleH [degrees]')

plt.subplot(313)
plt.plot(D[0:10000])
plt.title('LAS T3D_Distance over 10000 points')
plt.ylabel('Distance in meters')
plt.xlabel('Integer Point Record')

plt.show()